Optimal Multilevel Methods for H(grad), H(curl), and H(div) Systems on Graded and Unstructured Grids
نویسندگان
چکیده
We give an overview of multilevel methods, such as V-cycle multigrid and BPX preconditioner, for solving various partial differential equations (including H(grad), H(curl) and H(div) systems) on quasi-uniform meshes and extend them to graded meshes and completely unstructured grids. We first discuss the classical multigrid theory on the basis of the method of subspace correction of Xu and a key identity of Xu and Zikatanov. We next extend the classical multilevel methods in H(grad) to graded bisection grids upon employing the decomposition of bisection grids of Chen, Nochetto, and Xu. We finally discuss a class of multilevel preconditioners developed by Hiptmair and Xu for problems discretized on unstructured grids and extend them to H(curl) and H(div) systems over graded bisection grids.
منابع مشابه
Algebraic Multigrid for High-Order Hierarchical H(curl) Finite Elements
Classic multigrid methods are often not directly applicable to nonelliptic problems such as curl-type partial differential equations (PDEs). Curl-curl PDEs require specialized smoothers that are compatible with the gradient-like (near) null space. Moreover, recent developments have focused on replicating the grad-curl-div de Rham complex in a multilevel hierarchy through smoothed aggregation ba...
متن کاملFinite Element Methods for Maxwell Equations
1. SOBOLEV SPACES AND WEAK FORMULATIONS Let Ω be a bounded Lipschitz domain in R. We introduce the Sobolev spaces H(curl ; Ω) = {v ∈ L(Ω), curlv ∈ L(Ω)}, H(div; Ω) = {v ∈ L(Ω),div v ∈ L(Ω)} The vector fields (E,H) belong to H(curl ; Ω) while the flux (D,B) in H(div; Ω). We shall use the unified notation H(d; Ω) with d = grad , curl , or div. Note that H(grad ; Ω) is the familiar H(Ω) space. The...
متن کاملA nonconforming finite element method for a two-dimensional curl-curl and grad-div problem
Abstract. A numerical method for a two-dimensional curl-curl and grad-div problem is studied in this paper. It is based on a discretization using weakly continuous P1 vector fields and includes two consistency terms involving the jumps of the vector fields across element boundaries. Optimal convergence rates (up to an arbitrary positive ) in both the energy norm and the L2 norm are established ...
متن کاملLemme — — — - Weak vector and scalar potentials . Applications to Poincaré ’ s theorem and Korn ’ s inequality in Sobolev spaces with negative exponents .
In this paper, we present several results concerning vector potentials and scalar potentials with data in Sobolev spaces with negative exponents, in a not necessarily simply-connected, three-dimensional domain. We then apply these results to Poincaré’s theorem and to Korn’s inequality. 1 Weak versions of a classical theorem of Poincaré In this work, (the results of which were announced in [2]),...
متن کاملWavelet bases in H(div) and H(curl)
Some years ago, compactly supported divergence-free wavelets were constructed which also gave rise to a stable (biorthogonal) wavelet splitting of H(div; Ω). These bases have successfully been used both in the analysis and numerical treatment of the Stokes and Navier–Stokes equations. In this paper, we construct stable wavelet bases for the stream function spaces H(curl; Ω). Moreover, curl-free...
متن کامل